  do is = 0, ns - 1

     print("is : "+is) 

     if(sitename(is) .eq. "mauna_loa")
       lonx(is) = lonx(is) + 30
       latx(is) = latx(is) 
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "bermuda")
       lonx(is) = lonx(is) + 27
       latx(is) = latx(is) 
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "livermore")
       lonx(is) = lonx(is) - 7
       latx(is) = latx(is) 
       kset(is) = 0 
     end if
                
     if(sitename(is) .eq. "socorro")
       lonx(is) = lonx(is) - 4  
       latx(is) = latx(is) - 2 
       kset(is) = 0 
     end if
                
     if(sitename(is) .eq. "grifton0001")
       lonx(is) = lonx(is)   
       latx(is) = latx(is) - 5
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "para")
       lonx(is) = lonx(is)  
       latx(is) = latx(is) 
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "chacaltaya")
       lonx(is) = lonx(is) - 9
       latx(is) = latx(is) 
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "rio_de_janeiro")
       lonx(is) = lonx(is) - 14
       latx(is) = latx(is) 
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "pathankot")
       lonx(is) = lonx(is) - 8
       latx(is) = latx(is) 
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "bombay")
       lonx(is) = lonx(is) - 5
       latx(is) = latx(is) 
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "cape_point")
       lonx(is) = lonx(is) - 8
       latx(is) = latx(is) 
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "crozet")
       lonx(is) = lonx(is) - 3
       latx(is) = latx(is) 
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "mawson")
       lonx(is) = lonx(is) - 5
       latx(is) = latx(is) 
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "kerguelen")
       lonx(is) = lonx(is) + 8 
       latx(is) = latx(is) - 7
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "amsterdam_island")
       lonx(is) = lonx(is)  
       latx(is) = latx(is) + 7 
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "cape_grim")
       lonx(is) = lonx(is) - 9 
       latx(is) = latx(is) 
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "macquarie_island")
       lonx(is) = lonx(is) - 18
       latx(is) = latx(is)
       kset(is) = 0
     end if

     if(sitename(is) .eq. "dumont")
       lonx(is) = lonx(is) - 5
       latx(is) = latx(is)  
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "hongkong")
       lonx(is) = lonx(is)  
       latx(is) = latx(is) - 5 
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "beijing")
       lonx(is) = lonx(is) + 22
       latx(is) = latx(is) 
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "changchun")
       lonx(is) = lonx(is) + 15 
       latx(is) = latx(is) + 5 
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "huhehaote")
       lonx(is) = lonx(is) - 9 
       latx(is) = latx(is) + 4
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "xi_an")
       lonx(is) = lonx(is) - 1
       latx(is) = latx(is) + 1 
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "guiyang")
       lonx(is) = lonx(is) - 4
       latx(is) = latx(is)  
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "fuzhou")
       lonx(is) = lonx(is) + 23
       latx(is) = latx(is) - 4
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "gosan")
       lonx(is) = lonx(is) + 22
       latx(is) = latx(is) + 1
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "shanghai")
       lonx(is) = lonx(is) + 26 
       latx(is) = latx(is) - 4
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "moscow")
       lonx(is) = lonx(is) + 25
       latx(is) = latx(is) - 3 
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "zingst")
       lonx(is) = lonx(is) + 21
       latx(is) = latx(is) - 3 
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "vaasa")
       lonx(is) = lonx(is) - 2 
       latx(is) = latx(is)  
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "pallas")
       lonx(is) = lonx(is) - 1
       latx(is) = latx(is) + 3 
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "kevo")
       lonx(is) = lonx(is) + 17
       latx(is) = latx(is) + 5 
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "rovaniemi")
       lonx(is) = lonx(is) + 29
       latx(is) = latx(is) + 2 
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "nurmijarvi")
       lonx(is) = lonx(is) + 15
       latx(is) = latx(is) - 3.6
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "joensuu")
       lonx(is) = lonx(is) + 25
       latx(is) = latx(is)  
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "saclay")
       lonx(is) = lonx(is) - 2
       latx(is) = latx(is)  
      ;kset(is) = 0 
     end if

     if(sitename(is) .eq. "toulouse")
       lonx(is) = lonx(is) - 6 
       latx(is) = latx(is)  
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "mace_head") 
       lonx(is) = lonx(is) - 10
       latx(is) = latx(is) + 2 
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "cincinnati") 
       lonx(is) = lonx(is) - 6
       latx(is) = latx(is) + 5 
       kset(is) = 0 
     end if

     if(sitename(is) .eq. "hohenpeissenberg")
       lonx(is) = lonx(is)  
       latx(is) = latx(is)  
       ;prtname(is) = "2" 
     end if

     if(sitename(is) .eq. "heidelberg")
       lonx(is) = lonx(is)  
       latx(is) = latx(is)  
       ;prtname(is) = "3" 
     end if

     if(sitename(is) .eq. "schauinsland")
       lonx(is) = lonx(is)  
       latx(is) = latx(is)  
       ;prtname(is) = "4" 
     end if

     if(sitename(is) .eq. "jungfraujoch")
       lonx(is) = lonx(is)  
       latx(is) = latx(is)  
       ;prtname(is) = "5" 
     end if

     if(sitename(is) .eq. "freiburg")
       lonx(is) = lonx(is)  
       latx(is) = latx(is)  
       ;prtname(is) = "6" 
     end if

     if(sitename(is) .eq. "gif_sur_yvette")
       lonx(is) = lonx(is) - 13  
       latx(is) = latx(is)  
       kset(is) = 0 
       ;prtname(is) = "7" 
     end if 

     if(sitename(is) .eq. "lutjewad")
       lonx(is) = lonx(is)  
       latx(is) = latx(is)  
       ;prtname(is) = "8" 
     end if

     if(sitename(is) .eq. "paris")
       lonx(is) = lonx(is)  
       latx(is) = latx(is)  
       ;prtname(is) = "9" 
     end if

     if(sitename(is) .eq. "new_york")
       lonx(is) = lonx(is)  
       latx(is) = latx(is)  
       ;prtname(is) = "10" 
     end if

     if(sitename(is) .eq. "chester")
       lonx(is) = lonx(is)  
       latx(is) = latx(is)  
       ;prtname(is) = "11" 
     end if

     if(sitename(is) .eq. "washington_dc")
       lonx(is) = lonx(is)  
       latx(is) = latx(is)  
       ;prtname(is) = "12" 
     end if

     if(sitename(is) .eq. "sterling")
       lonx(is) = lonx(is)  
       latx(is) = latx(is)  
       ;prtname(is) = "13" 
     end if

     if(sitename(is) .eq. "nanjing")
       lonx(is) = lonx(is)  
       latx(is) = latx(is)  
       ;prtname(is) = "14" 
     end if

     if(sitename(is) .eq. "wuhan")
       lonx(is) = lonx(is)  
       latx(is) = latx(is)  
       ;prtname(is) = "15" 
     end if

     if(sitename(is) .eq. "milan")
       lonx(is) = lonx(is)  
       latx(is) = latx(is)  
       ;prtname(is) = "16" 
     end if

     if(sitename(is) .eq. "erba")
       lonx(is) = lonx(is)  
       latx(is) = latx(is)  
       ;prtname(is) = "17" 
     end if 

  end do











